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SUMMARY 


As part of an . investigation of flows around and forces on an oscil- 
lating finite wing in subsonic compressible flow, the exact double inte,- „ 
grals occurring in the theoretical formulation of the three-dimensional 
nonstationary aerodynamic forces have previously been reduced to single 
integrals over the range of either one of the two independent variables. 
The present report describes a method for solving the resultant three- 
dimensional problem based on a generalization of approximate methods for 
solving two-dimensional problems.' It is shown that the calculation of 
three-dimensional corrections to the two-dimensional theory involves 
only the solution of a one -dimensional integral equation for the spanwise 
variation of circulation, provided tabulated values for the kernel of the 
integral equation are available. 


INTRODUCTION 

The present report is a sequel to an earlier report (reference 1) 
on the problem of the oscillating lifting surface of finite span in sub- 
sonic compressible flow. In the earlier report the exact double -integral 
equation of the problem had been reduced to an approximate integral equa- 
tion in such a way that the double integrals are reduced to single inte- 
grals over the range of either one of the two independent variables. The 
present report describes a method of solution of the resultant three- 
dimensional problem. The proposed method is a generalization of any one 
of the known approximate methods for the solution of the two-dimensional 
problem. It is shown that, the calculation of three-dimensional correc- 
tions to the two-dimensional theory involves, within the framework of the 
proposed procedure, nothing more than the solution of a one -dimensional 
integral equation for the spanwise variation of circulation. The solu- 
tion of this one -dimensional integral equation will be no more difficult 
than the corresponding solution for the problem of incompressible flow, 
once the kernel of the integral equation has been tabulated for an appro- 
priate range of parameters. Determination of the numerical values’ of 
this kernel function is not a straightforward matter because of the 
particular analytical form in which this function appears. Applications 
of the present theory depend on the execution of this numerical task. 
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SYMBOLS 


x* 

y* 

G,K 

V* = b*V 


x m* 


M2) 


fi* 

T,* 

b* 

x,y,z 

k 

a 

g . 

X. 

V = 


= Tj/s \J 1 - 


1 - M 2 


Pa 


quantities defined as g*(x*,y*) = g(x,y); 
X*(x*,y*) = X(x>y) 

dimensionless coordinate ((x - x m *)/b*) 
dimensionless spanwise coordinate (y/s \[T 
dimensionless coordinate 
kernels of integral equation (l) 


quantity indicating amount of sweep 
quantity defined by equation (2) 

M 2 

local semichord divided by semichord b at midspan 

dimensionless coordinates 

reduced-frequency parameter 

circulation function 

function defined by equation (8) 

function defined by equation (6) 

variables of integration 

pressure amplitude at lifting surface 
modified potential function 
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ta = t(x,y,0) 


Po 

U 



density of stream flowing with velocity U 
main-stream velocity 
semichord at midspan 


H 

M 

a 

i 

p 

Q.,v± 

Sm 

S 

Ho ( 2 W 2) 

f m 

o,t 

Ap a 



p(2) 

B 

I,H,l2 


instantaneous shape of lifting surface 

Mach number of main stream 

arbitrary constant 

auxiliary variable of integration 

quantity defined by equation ( IT) 

functions defined by equations ( 19) and (20) 

function defined by equation (22) 

ratio of semispan to semichord at midspan 

Haiikel functions of second kind and. of zeroth and first 
order y respectively 

function defined by equation ( 23 ) 

auxiliary variables of integration 

component of pressure distribution p a representing 
solution of special two-dimensional problem 

component of pressure distribution ,p a representing 
solution of problem when three-dimensional effects are 
neglected 

quantity defined by equation (28) 
quantity defined by equation (a 4) 
functions defined by equations (A5) 
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JO,Y 0 

Fi,F 2 


arbitrary constant 

Bessel functions of zero order 

real and imaginary parts, respectively, of function F M 
defined by equations (A12) 


FORMULATION OF PROBLEM 


The problem which. is considered in what follows consists in the 
solution of an integral equation of the following form: 


g*(x*,y*) = * <f X*(|*,y*)G(x* - 6*) d£* + 


1-1 


iv*e _iV * Xm * ft*(y*) e" iv *^* G(x* - £*) d|* + 

1 


s -iv*( x*+xm*) 


1-1 


dft* 

dr,* 


K( y* - r,*) dq* 


( 1 ) 


The form of the functions G and K follows from a comparison of the 
foregoing equation with equation (88) of reference 1. Equation (l) is 
to be solved for X*, in terms of g* and ft*, and then use is to be 
made of the definition of ft* in terms of X*; namely. 


ft* = e iv 


X* dx* 


( 2 ) 


1-1 


Attention shall here be restricted to the problem of the rectangula 
lifting surface for which 


( 3 ) 


b* = 1 

v* = 

X* = X 

ft* = 

x m * = 0 . 

g* = 
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Equation (l) may then be rewritten in the form 

g(x,y*) = - J* X(|,y*)G(x - |) dl + 

n oo 

' ' ivfi(y*). / e -iv ^ G(x - £) d| + 


5 _ivx i r K(y* - q*) dr,* 

ar|* 

-i 


w 


n(y*) = e lV j X(x,y*) dx. 

-1 


( 5 ) 


It is now recalled that the function X is related to the pressure 
amplitude p a at the lifting surface by means of the equations 


and 


d^x,y,0) _ p dt a 
dx dx 



(6) 


( 7 ) 


and that the function g is determined in terms of the instantaneous 
shape H of the lifting surface by the relation 



Ue" 1 ^ 

Vl - M 2 



(8) 


Note also that 





( 9 ) 


V - [i = k 
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The proposed procedure for the solution of equations ( 4 ) and ( 5 ) 
depends on the introduction of the function p a instead of X into 

the integral equation, that is, on the reduction of this equation, 
when dfi/dq* = 0, to Possio’s form.- This step will be carried out in 
the next section. 


INTEGRAL EQUATION FOR PRESSURE DISTRIBUTION AT AIRFOIL 


Combination of equations (7). and ( 9 ) leads to the following 
alternate relation between and p a : 


ty a U,y*) = - ttu e " ivx j P a ( x '>y*) eikx ' dx ' 

J-l 


do) 


A corresponding expression for the circulation function Q is obtained 
as follows: 


iv 


ft = 2e lv I = 2e iV 


dx 


1-1 


= -^U / P a (x , ,y*)e lkx ' dx* 

"o u 

-1 


(ID 


From equations (6) and (10) it follows that 


X = - 


2b 

PnU 


i(k-v)x — . -ivx I — ikx' . . 

e . ' p a - ive p a e dx 1 

-1 


( 12 ) 


Equations (ll) and (12) are introduced into integral equation ( 4 ). 

There occurs one term which must be transformed by integration by parts, 
namely, the term 
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e" iv ^ G(x - |) p a e ikx ’ dx' | d| = 

l-l[_ J-i 

i pi 

e-ivl' G(x- - •!') . d|* I P a elkX ' dx r 

a J-l 

1 


1-1 


r* 

^i 

e“ iv i' G(x - I’ ) d£ ' 

p a (i^y*) elk| > 

L 



a 



y d | 


(13) 


The arbitrary constant cl may be set equal to infinity. If this is 
done then the integrated portion of equation (13) can be used to cancel 
the second term in integral equation (4). Equation (4) becomes 


g = 


2b 

PqU 


p a U*y*) 

1 


e i(k-v) I q( x - i) + 


iv 


e" iv| ’ G(x - I’) dl’e 1 * 1 


d| + 


pl 


-ivx 


d n 

dq* 


K( y* — T|*) dr,* 


(14) 


-l 


On the left-hand side of equation (l4) the value of g from equa- 
tion (8) is introduced. Then both sides of the resultant equation are 


multiplied by (poU^bJe 1 ^. This gives 

* If)- f 

2b \f 1 - \ 


1-1 


;ip(x-i) G(x - |) - 


ive 


ipx+ikl I e - iv i * G ( x _• |.) d |« 


'6 

■>1 


d| + 


PpU -ikx 
2b 6 


dft 

dr,* 


K( y* - . r, * ) dr,* 


(15) 


-1 
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Now set 


>x-£ 


G(x : I) = - I) - ive ik ( x ^) e iv C G( £) d£ (16) 


and 


Pi 


Po u 

2b 


ft = - 


p a ( x >y*) elkx = - p 


(17) 


This leads to the following final form of the integral equation of the 
problem ■' 

= <f Pa( &>y*)G(x - i) - 


2b 


f lkB ♦ jS^ - 
rr# v 3x ' 


75 


pi 


-ikx 


-1 


dP 

dri* 


K(y* - q*) d^* 


(18) 


When dP/drj* = 0 this integral equation must reduce to the corresponding 
equation of the two-dimensional theory; that is, the kernel G must he 
reducible to the corresponding kernel in Possio’s formulation of the 
problem. It is not necessary for the present purposes to verify explic- 
itly that this requirement is satisfied. 


It is useful to introduce as further abbreviations 

nl 


*1 


-1 


dP 

dq* 


K(y* - q*) dq* 


(19) 


w i = 


Po^ 


2b\/l - M 2 


i« + p 

OX i 


(20) 
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Finally, before proceeding, the analytical expression for the 
kernel function K is listed explicitly. 


K(y* - r,*) = % 

kns VI - M2 


f l - M2 


(y* _ n*) 


where 


S M (z) * 

2 V 1 _ M 2 z 


»-M|z| 


Hq ( 2) (UI) <*£ 


F m (z) 


1 - M c 


FmU) 


LLL e _ia 

Z 


ae -iM N/a2 + ^2 ■ i 


- i 2 Wa2 + £ 


t 2 . 1 2 


+ iM ] d£ + 


-00 \ U “*°° 


,-iM VT^+f 


] t 2 + (2. 


d£ dT + 


iM \/t 2 +z 2 


1 2 I e 

T 2 + Z 2 


s/t^ + Z 2 


+ iM \ dT dcr 


The difficulty of the problem of evaluating the function K lies 
in the form of the function F»*. In an appendix to this report there 
(jre incorporated some preliminary considerations by Z. Kopax cuu^cming 
evaluation of this function Fpp 


SOLUTION OF INTEGRAL EQUATION FOR PRESSURE DISTRIBUTION p a 


Equation (l8) is written in the form 


Vi(x,y*) + e“ ikx Q(y*) =. <f P a ( Z>y*)G 4£ . 
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and it is seen that the expression for p a may he considered to he 
composed of two parts, as follows: 


P a (x,y*) = pj 2 )(x,y*) - Q Ap a (x) _ (25) 


_ ( 2 ) 

The term p a represents the solution of the problem when three- 
dimensional effects are neglected. The term Ap a represents the solu- 
tion of a special two-dimensional problem, namely, that problem for 


which w i (x,y*') = -e"- 5 -^*. This result holds regardless of the manner in 
which the solution of the two-dimensional problem is approached, and 
for the present purposes p a (2) and Ap a may be considered as known. 


It remains to determine the values of Q. • In order to do this 
equation (25) may be combined with equation (17). This gives' 


n 1 


P(y*) = 


p a ( 2 >e ikx dx 


-1 



(26) 


Use is now made of the definition of Q by means of equation (2 6 ), in 
order to write » 


P(y*) = 


Pi 


p a (2)e ikx dx - 


-1 


PI 


Ap e ikx dx 

cl 


dP 

dr,* 


K dr,* (27) 


-1 


1-1 


In equation (27) the abbreviations 

il 

p(2)(y*) _ I p .( 2) ( x ,y*) e ikx ^ 


(28) 


1-1 


il 

pi(k,M) =• Ap a e ikx dx 

J-l 


(29) 


are used in order to obtain the basic integral equation of the problem 
in the form 

^1 

P(y*) + n(k,M) <j ^ 5 - K(y* - r,*) dq* = P ( 2 ) (y*) (30) 

J-l 
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Equation (30) is the appropriate generalization of the previously- 
obtained integral equation for the circulation function fi in incom- 
pressible flow (reference 2). It is seen that determination of the 
terms p( 2 ) and p, as defined by equations (28) and (29), requires the 
solution of a two-dimensional problem in accordance with equations (24) 
and (25). In view of equations (19) and (30) the factor Q in the 
three-dimensional correction term in equation (25) may also be written in 
the form 



(3D 


What is left to be done, in order to obtain three-dimensional 
corrections of the kind here contemplated, is to solve integral equa- 
tion (30), and there the main difficulty is the tabulation of the func- 
tion K. Once. K has been tabulated the problem of subsonic compress- - 
ible flow is solved to the same extent as the problem of incompressible 
flow has been solved in reference 2. The solution of specific problems 
may then be obtained in the same manner as was done in reference 3 Tor 
problems of incompressible flow. 


Massachusetts Institute of Technology 
Cambridge, Mass., May 5, 19^9 
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APPENDIX 

SOME CONSIDERATIONS ON EVALUATION OF FUNCTION F M 
By Z. Kopal 


Given 


F m (x) 


= M 


poo 


-i o 


^ ^ ae -iM \/<J 2 + I 2 


+ t 2 \ n/ct 2 + 5 2 


+ iM\ d| + 


^-imn/^ 2 \ 

- — ' f - + iM ] dT ■+ 

T 2 + X 2 \ \/t 2 + X 2 / 






dr d| 


F —oo vy -oo 


da 


(Al) 


Since 


1- x 


-iMs/a^S 2 . A 

Wns * ”) 


«-s 


,-lM 




l| x | n/o 2 * 


d| 


+ 5 £ 


|x|e -iM\/^ , ! 

t2 + * 2 V \FTIs 


IM i dT = 


|x| d_ 

x dx' 


poo 


x/t^h 


,-iM NJt^+x 


7 


dT 


T 2 + X 2 


(A2) 


and since, moreover, the order of integration in the last integral 
(factored by M 2 ) is interchangeable, the foregoing function F M (x) can 
also be written as 


Fj,j(x) = - e - -*- 0 B(a,x;M) da 


(A3) 
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where 


noo 


B = l(x,o2)- + I(ff,x2) - M 2 


I(ff,t 2 )-d«y 


(AU) 


In equation (A^) the abbreviations 
l(a,a 2 ) s 


; -iM \/a 2 +y 2 

\}W 


dy = I^a-jCt?) - il 2 (a,a 2 ) 


+ a* 


/ o\ I cos M V a 2 + y 2 , 
Il(a,a 2 ) = f= / dy 


fa 


r>°° 


I 2 (a,a 2 ) = 


\/a£ + y2 


sin Mv/a 2 + y 2 
\|a2 + y 2 


dy 


(A5) 


are used and (since only positive values of x are of interest here) 
|x|/x is set equal to 1. 

The whole difficulty of the problem centers around the evaluation 
of the I integrals when M / 0. If a = 0, 


0 


£ ° s , M\/a 2 + £ dy = g y (jfa) 
\la? + y 2 2 


SlnMN/ ? 2 . l y2 dy = f J 0 (Ma) 


(A6) 




It can, f urthermore , be proved (although the proof is somewhat involved) 
that the I functions, regarded as functions of M, satisfy the following 
differential equations: 
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M 2 11 + M + aPvP-I-, = Ma sin M\/a 2 + a 2 

aw 2 * 1 


^ L?- + M + a^M^l2 = Ma cos M s/a^ + a? 

Sm 2 3m 


(A7) 


which, incidentally, can be written as a single equation of the form 


M 2 + M & + a 2 M 2 I = -iMae 
&M 2 ® 


iMx/a^+a' 


2 ' 


(A8) 


If a = 0, the homogenous part of the foregoing equations is easily 
recognized as Bessel’s equation of zero order and their solutions are 


Il = | Y Q (Ma) 


I 2 = | J 0 (Ma) 


I = II - il2 ; = -#- Ho (l) (Ma) 


(A9) 


where IL^ 1 ) denotes the respective Hankel function.. 

When the two particular solutions of the homogeneous equation 
(for a =0) are known, the inhomogeneous equations admit of a solution 


PM 


I-L = A y 0 (Mcl) + aJ 0 (Ma) 


| 2 Y 0 Ua) sin £ 


\[a? + 


<4 


nM 


aY 0 (Ma) 


| 2 J 0 ( |a) sin | \J a 2 + a 2 d£ 


(A10) 


and 
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>M 


I 2 = | J 0 (Ma) + aY 0 (Ma) I | 2 J 0 (£a) cos | \J a 2 + a 2 d| - 




ajQ(Ma) l 2 Yo(la>) cos 5 \/a 2 + a 2 d£ 


Since both transcendent functions occurring behind the integral 
sign of these integrals can be expanded in rapidly converging power 
series of i and subsequently evaluated, these integrals can hence- 
forward be regarded as known. This leaves one more quadrature to be 
performed in order to evaluate the function F^(x) - the real and 
imaginary parts of which can easily be written as 


F m (x) = F x - iF 2 


where 


** 

II 

"'l 00 

1; Jl(x^ a 2 ) + I]_( a,x 2 ) M 2 1 

Y CO 

Il(a,! 2 ) cL6 

J 

0 

J 

X 


cos a da - 


2- I 2 (x,a 2 ) + — I 2 (a,x 2 ) - M 2 I 2 (a,| 2 ) d| 
da ox | 

X 


sin a da 


F 2 = 


B I 1 (x,a 2 ) + |- I x ( cr,x 2 ) 


do 


- Vf 


Il(a,| 2 ) d| 


sin a da + 


p 00 


o L 


I 2 (x,a 2 ) + I 2 (a,x 2 ) 


dx 


- M 2 I 2 ( a,|2) d| 


cos a da 


(All) 


>(A12) 
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